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ABSTRACT 

By the time the expanding envelope of a Type II supernova becomes transparent in the optical continuum, 
most of the gamma-ray luminosity produced by radioactive Fe/Co/Ni clumps propagates into the hydrogen/ 
helium envelope and is deposited there, if at all The resulting fast electrons excite He i and H i, the two- 
photon continua of which are the dominant internal sources of ultraviolet radiation. The UV radiation is 
blocked by scattering in thousands of resonance lines of metals and converted by fluorescence into optical and 
infrared emission lines that escape freely. We describe results of Monte Carlo calculations that simulate 
non-LTE scattering and fluorescence in more than five million allowed lines of Ca, Sc, Ti, V, Cr, Mn, Fe, Co, 
and Ni. For a model approximating conditions in the envelope of SN 1987A, the calculated emergent spec- 
trum resembles the observed one. For the first 2 yr after explosion, the ultraviolet radiation (A < 3000 A) is 
largely blocked and converted into a quasi continuum of many thousands of weak optical and infrared emis- 
sion lines and some prominent emission features, such as the Ca ii AA8600 triplet. Later, as the envelope cools 
and expands, it becomes more transparent, and an increasing fraction of the luminosity emerges in the UV 
band. 

Subject headings: line: formation — radiative transfer — stars: atmospheres — supernovae: general — 
supernovae: individual (SN 1987A) 


1. INTRODUCTION 

A major challenge in interpreting the spectra of many astro- 
physical systems is to understand the blocking of ultraviolet 
(UV) radiation by the resonance lines of metals. For example, 
UV line blocking is known to play a major role in spectrum 
formation in the atmospheres and winds of hot stars (Abbott & 
Lucy 1985; Kudritzki 1989), active galactic nuclei (Netzer & 
Wills 1983), classical novae (Hauschildt et al. 1992), and super- 
novae (Karp et al. 1977; Lucy 1987; Wagoner, Perez, & Vasu 
1991; Hauschildt & Ensman 1994). The success of the Interna- 
tional Ultraviolet Explorer and the Hubble Space Telescope in 
routinely observing the UV spectra of such systems has 
increased the stakes in understanding this UV line blocking. 

We are especially motivated to address this problem for 
supernova envelopes for three reasons. First, SN 1987 A is the 
first supernova for which excellent UV, optical, and infrared 
spectra are available well into the nebular phase, defined as the 
time when the photosphere has apparently vanished and the 
envelope has become transparent in the optical continuum. SN 
1987A entered its nebular phase at t « 4 months. Thereafter, 
the optical spectrum was dominated by emission lines and the 
fraction of the radiation that escaped as UV radiation 
increased (McCray 1993). Indeed, SN 1987A has provided the 
first detailed view of the evolution of the UV spectrum of a 
supernova, from outburst to the present time (Pun et al. 1995). 

Second, we have been engaged in a systematic effort to 
model the evolution of the spectrum of SN 1987A. In doing so, 
we have realized that a UV continuum field within the 
envelope is required to account for the observed light curves of 
some of the most prominent optical and infrared emission lines 
such as the hydrogen recombination lines (Xu et al. 1992), the 
Ca ii infrared triplet (Li & McCray 1993), and the Fe n A7157 
lines (Li, McCray, & Sunyaev 1993). From such lines we can 
infer the intensity of this UV radiation field, even during the 


first 2 yr after explosion, when line blocking prevents most UV 
radiation from emerging. 

Third, the outstanding unsolved problem in understanding 
the nebular spectrum of SN 1987 A is to account for the origin 
of the optical continuum that underlies the emission-line spec- 
trum. At the temperatures inferred from the optical and infra- 
red emission lines, the envelope is too cool and neutral to emit 
this continuum by any known thermal processes, such as 
bremsstrahlung or recombination (Wooden et al. 1993). Here 
we show that this continuum is actually many thousands of 
overlapping fluorescence lines resulting from the resonant 
absorption of UV continuum radiation in the envelope. Thus, 
understanding this fluorescence is the last major link in the 
logical chain that explains how the luminosity of gamma rays 
is converted, first to the energy of fast electrons, then to atomic 
excitation and heating, then (partly) to the UV continuum, and 
finally to the UV, optical, and infrared spectrum that is 
observed. In § 2 we describe the basic physics controlling the 
resonant scattering and fluorescence of UV photons propagat- 
ing through an expanding medium and discuss previous work 
on this problem. In § 3 we describe our Monte Carlo technique 
for solving this problem. (In Appendix A we test this technique 
with an idealized problem that can be solved analytically.) In 
§ 4 we describe a more realistic model for a supernova 
envelope including a source of UV continuum embedded in a 
hydrogen envelope with cosmic (Large Magellanic Cloud) 
abundances. As we show in § 5, the UV and optical light curves 
and spectra that emerge from such a model resemble those 
observed from SN 1987 A. 

We emphasize, however, that the present paper is only a step 
on the way to a complete theory for the emergent optical/UV 
spectrum of SN 1987A (and other supernovae). Before attempt- 
ing to model the observations in detail, we must include a 
number of physical effects not considered here. In § 6 we sum- 
marize the progress to date and outline the work to be done. 
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2. LINE TRANSFER IN A SUPERNOVA ENVELOPE 

The problem of line transfer in a freely expanding supernova 
envelope can be described quite accurately with the Sobolev 
approximation (Sobolev 1957, 1960; Castor 1970; Mihalas 
1978). Here we briefly summarize the method. For further 
details, see Xu et al. (1992), Li & McCray (1992, 1993, 1995), 
and Li et al. (1993). In a medium that is expanding supersoni- 
cally (with velocity t> eKp ), any photon is resonant to a particular 
atomic transition only in a layer of thickness 


Sr 



( 1 ) 


where v th % (kT/m) 112 . For example, for Fe ii lines at 7000 K, 
t? lh — 1 km s' 1 . Within this resonant layer the photon may be 
scattered by the resonant line many times before it escapes 
from the layer. To account for resonant trapping in the 
Sobolev approximation, one should replace the normal spon- 
taneous transition rate A ul (Einstein A) by the effective tran- 
sition rate 

= A ul P C J t J . (2) 

where the subscripts u and l represent the upper and lower 
level, respectively. The Sobolev escape probability is given by 


Pc sc( T J = 


1 - exp (-tJ 


( 3 ) 


and the Sobolev optical depth of the line by 
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where f lu is the absorption oscillator strength, A 0 is the line 
wavelength, the gs and N's are the statistical weights and 
population densities of the relevant levels, respectively, t is the 
time since explosion, and the term in parentheses accounts for 
induced emission. 

As the photon leaves the resonant region, its wavelength is 
redshifted relative to the comoving gas according to Hubble's 
law S?./a 0 = H 0 Sr/c, where the “ Hubble constant ” is given by 


// _ -L ^ p. _ t 
H 0 - ^ -t 


( 5 ) 


Accordingly, the photon may come into resonance with 
another atomic transition at a longer wavelength. Therefore, 
for a photon of given wavelength created at any location in the 
envelope, the resonance transitions to the red of that photon 
wavelength define a sequence of concentric spherical 
“resonance shells." Figure 1 illustrates the radiative transfer. If 
the photon is scattered at a given resonance shell, a new 
photon is emitted with random direction at that point. Once 
again, the resonance transitions to the red of the new photon 
define a new set of resonance shells, and the process repeats 
until the photon escapes the emitting region. 

In fact, it is likely that the atom excited by resonance absorp- 
tion will not reemit the original photon, but instead will 
cascade down by emitting several photons with longer wave- 
lengths (fluorescence). The atom may also be excited or de- 
excited by thermal electron collisions. The branching ratio for 
a particular radiative transition from upper state u to state i is 



Fig. 1. — Propagation of a photon that experiences multiple line scattering 
in an expanding medium. 


given by 

A ef{ 

Pui = yj< + wii ’ (6) 

where C UJ {T ) is the collisional excitation or de-excitation rate 
coefficient, and we have used the notation Aff = 0 for j < i for 
convenience. 

During the nebular phase, the ionization and level popu- 
lations of the metals responsible for most of the resonance line 
opacity are far from local thermodynamic equilibrium (LTE). 
Line blocking due to resonance transitions of metals renders 
the envelope largely opaque at UV wavelengths. As we shall 
see, fluorescence plays a dominant role in the escape of UV 
radiation from the envelope. UV photons, produced by the 
degradation of gamma rays, are largely converted into visible 
and infrared photons, to which the envelope is largely trans- 
parent. 


3. MONTE CARLO METHOD 

We model the radiative transfer and fluorescence of UV 
radiation through a supernova envelope with a Monte Carlo 
simulation. The details are described below and illustrated by 
the sketch in Figure 1. 

1. Choose the wavelength of the initial parent photon 
according to the prescription 

^0 = ^min “L 0i(-^max ^min) ’ 

where 0j and all the other 0/s in this section are random 
numers uniformly distributed on the interval 0 < 0, < 1. The 
source spectrum has lower bound A min = 912 A, because UV 
source photons with shorter wavelength will be photoabsorbed 
by H i in the envelope, which is nearly neutral (Xu et al. 1992). 
We choose an upper bound at / max = 5912 A because the 
envelope is nearly transparent to all radiation at longer wave- 
lengths. To account for such photons in the source spectrum, 
we can simply add them to the emergent spectrum. 

2. Give the initial photon a weight, w 0 , proportional to the 
emission probability of the source mechanism. For He i two- 
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photon emission, the w 0 can be calculated from the spectrum 
shape given by equation (14). 

3. Choose the radius at which the initial parent photon is 
emitted. In this paper we consider two models. In the first, all 
photons come from a central source. In the second, the source 
is uniformly distributed within a sphere, and so we choose a 
fractional radius according to 

do = ( 4 > 2) ,/3 . ( 8 ) 

where (f> 2 is another random number. 

4. Choose the angle 6 between the propagation direction of 
the initial photon and the radius vector to the point where it 
was produced, according to cos 6 = /i = 1 — 2</> 3 . 

5. Now calculate the distance from the point of origin to the 
surface of the envelope, given by 

r, = v /l - (1 - fi 2 )do - fid 0 . (9) 

A photon with original wavelength A 0 that propagates to the 
surface without scattering will be redshifted by A2 0 = 
X 0 r s v ma Jc relative to the comoving gas at the surface, where 
v max is the radial expansion velocity at the surface of the 
envelope. 

6. Calculate the net Sobolev optical depth, t s , along the ray 
r s to the surface, given by 

= It,. (10) 

where the resonance lines are listed sequentially in order of 
increasing wavelength, the index k stands for the first line to 
the red of / 0 , and the index s stands for the first resonance line 
(if any) to the red of 2 0 -I- AA 0 . Add a photon with wavelength 
A 0 and weight Wj = w 0 exp ( — i s ) to the appropriate wave- 
length bin of the emergent spectrum. 

7. For future reference, accumulate the weights, vvj, of each 
emergent photon,;, in two other arrays, one binned according 
to the number, n, of resonant scatterings that have taken place 
before the photon emerges, and the other binned according to 
the radius, r p of the last scattering event. 

8. If t s is negligible (say, t s < 0.01), return to step (1). Other- 
wise, give the remaining photon fraction a weight w 2 = w 0 
— w l — w 0 [l — exp ( — r s )]. This remaining fraction must be 
absorbed on one of the resonance shells along r s . To decide 
which one, choose an optical depth x x from the prescription 

1 , = -ln[l-*-*)*J- (11) 

(Note that eq. [1 1] ensures that 0 < x x < x s .) The photon will 
be absorbed by the first transition m for which 

m 

t,< It,. (12) 

i = k 

9. A resonance absorption creates an excited atom at a 
radius d x = (dl + r\ + 2 d 0 r x /*) 1/2 , where r x = (A m - A 0 )c/ 
(Aq v m ) (see Fig. 1). The atom (in upper level u) may decay to the 
initial state, or it may cascade through other intermediate 
levels. To decide, randomly choose a line j from among all m 
possible downward transitions according to 

" ’l < I A, , (13) 

i = 1 i = 1 i = 1 

where the fi ui 9 s are branching ratios given by equation (6). If the 
selected branch leaves the atom in an excited state, repeat the 


procedure to select subsequent branches until the atom relaxes 
to a metastable state. 

10. The cascade chosen in step (9) may create N new 
photons at radius d v Since it is inconvenient for a Monte 
Carlo code to follow all such photons, choose one of these 
photons at random and give it a weight w 3 = Nw 2 . This pro- 
cedure will not conserve energy in individual transitions, but it 
will do so in the aggregate with a sufficient number of trials. 

11. With this new photon and location, return to step (4) 
and repeat the process until no further resonant absorptions 
occur. Note that step (6) allows a fraction of the parent photon 
to escape with each resonance absorption. This procedure, 
which costs little extra computer time, vastly improves the 
statistics of the highly absorbed UV part of the emergent 
photon spectrum. 

Typically, we find that 20,000 initial continuum photons are 
sufficient to provide an emergent spectrum with reliable 
counting statistics throughout. 

4. MODEL SUPERNOVA ENVELOPE 

The envelope of SN 1987A consists of primordial gases 
mixed inhomogeneously with fresh nucleosynthesis products. 
The nucleosynthesis products, especially the Fe/Co/Ni, cer- 
tainly play a role in the scattering of UV radiation and the 
formation of the emergent spectrum. However, in this work we 
shall attempt only to illustrate the transfer of UV radiation 
through an idealized model envelope, having uniform density 
and temperature and composed of hydrogen and helium with 
cosmic abundances of metals. Although this choice is made for 
simplicity, it may be a reasonable first approximation to the 
actual radiative transfer in the envelope of SN 1987 A. As Li et 
al. (1993) have shown, after about 200 days, more than 90% of 
the gamma-ray luminosity absorbed by the supernova 
envelope is deposited in the hydrogen and helium gases that 
comprise most of mass of the envelope. 

A much smaller fraction of the gamma-ray luminosity 
emerges as emission lines of hydrogen and helium, however. 
Where does the rest go? The answer is that some 25%-35% of 
the absorbed gamma-ray luminosity heats the gas, and most of 
the remaining fraction ionizes and excites He i and H i atoms 
(Xu & McCray 1991; Kozma & Fransson 1992). These atoms 
are largely prevented from decaying to the ground state by 
allowed transitions, owing to the very high optical depth of the 
resonance lines (Xu et al. 1992). Instead, they often decay by 
emitting two-photon continuum radiation. Indeed, Xu et al. 
(1992), Li & McCray (1992), and Li et al. (1993) inferred that a 
strong UV field must be present in the envelope of SN 1987 A 
to account for the light curves of some of the stronger optical 
and infrared emission lines. Most of these UV continuum 
photons do not escape the envelope, however. As we shall 
show, they are absorbed by resonance scattering by trace ele- 
ments and converted by fluorescence to optical radiation. 

In this paper, we illustrate line blocking and fluorescence of 
the He i two-photon continuum, the spectrum of which can be 
fitted by the expression (Drake, Victor, & Dalgarno 1969) 

'VO’He) = -VoU - I 1 - 23 'hJ 2,5 ) S~ 1 for 0 < y He < 1 , 

(14) 

where y He = v/v He2y is the fraction of the energy carried by one 
of the two photons, v Hc2y is the frequency difference between 
the He i 2s l S and He i Is l S states, and A Hc0 — 145 s _1 . 
Because H i strongly absorbs photons with wavelength less 
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than 912 A, we include only that part of the He i two-photon 
continuum with greater wavelength as a UV source in our 
model. 

We consider two models of the spatial distribution of the 
UV source: a uniformly distributed source and a central point 
source. Specifically, we take a hydrogen envelope with 
M(H) = 8 Af © (Li & McCray 1995) and maximum expansion 
velocity v max = 3000 km s" 1 . The envelope contains trace ele- 
ments with abundances 0.5 times the solar values (Lang 1980), 
typical for LMC planetary nebulae (Dufour 1984; Dopita & 
Meatherington 1991). We include in our code radiative tran- 
sitions due to nine elements (Ca, Sc, Ti, V, Cr, Mn, Fe, Co, and 
Ni), all of which we assume to be singly ionized owing to 
exothermic charge transfer with ionized hydrogen. The reso- 
nance lines of these ions should account for most of the UV 
opacity and fluorescence. (Other abundant ions, such as H i, 
He i, N i, N ii, C ii, O i, and O ii, have relatively few UV 
resonance lines.) 

We use a database of wavelengths and oscillator strengths 
for all allowed transitions of the relevant ions kindly provided 
by Kurucz (1994). The database includes a total of 5,223,718 
lines: 246 levels and 923 lines for Ca ii, 781 levels and 49,667 
lines for Sc ii, 3894 levels and 254,515 lines for Ti ii, 11,986 
levels and 836,510 lines for V n, 11,183 levels and 1,151,521 
lines for Cr ii, 13,436 levels and 687,800 lines for Mn ii, 10,921 


levels and 1,059,533 lines for Fe ii, 11,353 levels and 877,576 
lines for Co ii, and 9039 levels and 305,673 lines for Ni ii. 

For these allowed transitions the electron number density, 
/V e , in the envelope of SN 1987A is always less than the critical 
number density, N CT = Al\ { /C ul . Therefore, we can resonably 
neglect transitions caused by thermal electrons in our calcu- 
lations of UV radiation transfer and fluorescence. 

When an ion such as Fe ii is excited by UV resonance 
absorption, it will decay by radiative cascade until it reaches a 
metastable level, which must decay by emission of a forbidden 
line or by electron impact de-excitation. The Grotrian dia- 
grams for the Fe ii quartet states (Fig. 2) illustrate the great 
number of levels that may be excited by resonance transitions 
from the metastable levels, as shown in the lower panels. 
(Grotrian diagrams for the doublet and sextet states of Fe ii 
look similar.) In the Kurucz database, 39,720 (about 4%) of the 
allowed transitions of Fe n are coupled to metastable levels. 
These are the only lines that can contribute to the UV opacity, 
because only the lower lying metastable levels have significant 
populations in the relatively cool supernova envelope. The 
subordinate allowed transitions among nonmetastable levels 
play a crucial role in the fluorescence, however. 

To calculate the UV opacity of the supernova envelope, we 
need to know the populations of ions in the metastable levels. 
As a first approximation, we might assume that these levels are 
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Fig. 2.— Grotrian diagram of Fe n quartet states. Metastable and nonmetastable, odd and even states are shown in different windows. 
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populated in LTE, with excitation temperature equal to the 
electron kinetic temperature. We do not take this assumption 
for granted, however. For the case of Fe ii, which accounts for 
much of the UV opacity, we calculated non- LTE populations 
by solving coupled rate equations for entry and exit into each 
metastable level, including not only electron impact excitation 
and de-excitation, but also radiative decays by emission of 
forbidden lines, absorption of resonance lines, and radiative 
cascades into these levels following resonance absorption. 
Among the 10,921 energy levels of Fe n in the Kurucz data- 
base, the lowest 63 are metastable. There are 68 more metasta- 
ble levels at higher energy, but they do not contribute 
significantly to UV opacity and can be neglected. 

By including the entry and exit rates due to allowed tran- 
sitions, we make the problem nonlinear. The populations of 
the metastable states now depend on the absorption and fluo- 
rescence of the UV photons, which in turn depend on the 
populations. The stationary populations of the lowest 63 meta- 
stable states obey the equation 

N ^L (NeCji + Af) + G,j = £ NlN t C tJ + A?/) + Sj , (15) 

where Gj is the exit rate from level j due to resonance absorp- 
tion, and Sj is the entry rate to level j due to radiative decays 
from higher levels (metastable or not). We have the constraint 

I NjGj-ZSj. (16) 

j j 

We adopt the following procedure to solve the populations 
of the metastable states: (1) make a first guess at the popu- 
lations by solving equation (15) with Gj = Sj = 0 (see Li et al. 
1993 for details); (2) using these populations, simulate the radi- 
ative transfer and fluorescence of the UV radiation field with 
the Monte Carlo code; (3) use the results of step (2) to calculate 
the S/s and Gj s; (4) calculate new populations of metastable 
states including the results of step (3) in equation (15); (5) 
repeat the Monte Carlo simulation with the new rates. 

One can see that the nonlinearity might introduce serious 
computational problems, and there is no guarantee that the 
procedure will converge. However, we find in practice that one 
iteration is enough. The spectrum that results from step (5) is 
not discernably different from the spectrum that results from 
step (2). For example, we find that for Fe n at 300 days and any 
temperature under consideration, electron impact collisions 
dominate the metastable level populations and so LTE is a 
very good approximation. At 800 days and 1000 K, the first 13 
levels of Fe ii are still approximately in LTE, but the upper 50 
levels are substantially overpopulated by UV pumping. 
However, even with these increased level populations, almost 
all resonance lines originating from these upper levels have 
x lu <$ 1. In other words, the opacity of the envelope is domi- 
nated by resonance transitions originating from the lowest 13 
levels. 

5. THE EMERGENT SPECTRUM 

Figures 3 a and 3 b show the calculated escape probabilities 
for our model of the envelope of SN 1987A at 300 days, with 
assumed temperatures of 6000 and 2000 K, respectively. Each 
resonance line is plotted with velocity resolution Av = 3000 
km s' l . The solid lines show the escape probabilities due to all 
elements in the database, and the dotted lines show the escape 
probabilities due to Fe ii alone. We see that Fe ii accounts for 




Fe ii only. 

much of the UV opacity. Notable exceptions are the bands 
3000 ^ X ;$ 3700 A, and 4100 < X < 4400 A, where most of the 
opacity is due to Ti n, and the strong absorption features due 
to Ca ii /U3934, 3969 (the H, K lines) and the Ca ii >U8600 
infrared triplet. 

Note that the envelope is mostly transparent for optical and 
infrared wavelengths but largely opaque at UV wavelengths. 
Note also the strong temperature dependence of the opacity. 
For T = 6000 K the envelope is almost completely opaque for 
X <> 3700 A. But, at T = 2000 K transparent “ windows ” 
appear in the UV spectrum. This occurs because the resonance 
lines in these windows couple to metastable states at higher 
energies. Since the populations of these states decrease rapidly 
with decreasing temperature, so do the Sobolev optical depths 
of the resonance lines in these windows. 

An isolated example of this effect is the Ca ii /U8600 infrared 
triplet, which couples the 3 d 2 D metastable state to a 4p 2 P 
state at higher energy (Li & McCray 1993). Note that the 
optical depth of this triplet is much greater at T = 6000 K (Fig. 
3 a) than at T = 2000 K (Fig. 3b). One might hope to use the 
observed depths of such absorption lines as temperature diag- 
nostics, but the interpretation is not so simple, owing to fluo- 
rescence. For example, in the spectrum of SN 1987 A, Ca ii 
/U8600 is seen in emission , not absorption. That occurs because 
its upper level is pumped by absorption of U V radiation in the 
H, K lines (Li & McCray 1993). 

Of course, the UV optical depths would also decrease with 
time owing to the expansion of the envelope (if temperature 
remained constant, we would have t oc t~ 2 ); but this effect is 
not nearly as dramatic as the effect of decreasing temperature. 

Figure 4 shows the results of a typical Monte Carlo simula- 
tion. For this case we assumed a uniformly distributed He i 
two-photon continuum source, shown as the thin solid curve 
(note that it is actually curved). We also assumed that the He i 
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Fig. 4. — The emergent spectrum for a uniformly distributed He t two-photon continuum source (UHeS) at 300 days and 6000 K 


two-photon source had the intensity calculated by Li & 
McCray (1995). The envelope has atomic density appropriate 
for t — 300 days and an assumed temperature T = 6000 K, 
close to the value inferred by Li & McCray (1993) from the 
Ca ii infrared emission lines. The emergent spectrum is shown 
with bin width A),/ a = v m Jc = 1/100, approximately the line 
widths that would be observed. The figure shows that the 
source spectrum is strongly absorbed for A < 2700 A, but much 
of the absorbed luminosity reappears as many thousands of 
fluorescence lines, mainly in the optical and near-UV bands 
where the lines overlap in a quasi continuum, shown by the 
solid spectrum. 

In addition to the fluorescence spectrum, a hot hydrogen 
envelope will also emit line radiation due to thermal electron 
impact excitation. To estimate the relative importance of such 
thermal line emission, we calculated the intensity of forbidden 
lines from the metastable levels of Fe n. The dashed spectrum 
in Figure 4 includes these forbidden lines. We see that at 
T = 6000 K the optical spectrum is dominated by thermal 
emission, not fluorescence. (We have not calculated the 
thermal emission lines from other ions in the envelope.) 

Figure 5 illustrates how the emergent spectrum changes as 
the supernova envelope expands and cools. Figure 5 represents 
the same model as Figure 4, except that the time is now 800 
days and the temperature has dropped to 2000 K. Accordingly, 
the density of the envelope has decreased by a factor 0.053 
owing to free expansion. The intensity of the He i two-photon 
continuum has also decreased by a factor ^ 10 2 owing to the 


decreasing illumination of the envelope by gamma rays. The 
emergent spectrum of Figure 5 diffes from that of Figure 4 in 
two major ways. First, a much greater fraction of the UV 
radiation escapes from the supernova envelope. Second, the 
optical spectrum is dominated by fluorescence. The envelope is 
now too cool to emit optical lines by thermal processes. Note 
however, that thermal electrons can still excite some prominent 
infrared emission lines. 

The solid vertical marks at the tops of Figures 4 and 5 
indicate fluorescence features that stand out above the quasi 
continuum. The marks in the top row of each figure indicate 
features dominated by Fe n, and those in the second row indi- 
cate fluorescence features dominated by other elements. Tables 
1 and 2 list these features for day 300 (Fig. 4) and 800 (Fig. 5), 
respectively. Wavelengths are listed with resolution A A/ A — 1/ 
100, and the percentage contributions (if greater than 10%) of 
various elements are given in parentheses. The remaining frac- 
tions may be attributed to weak lines and the original He i 
two-photon continuum. 

Note that the Ca ii AA8600 infrared triplet appears in emis- 
sion in both spectra. In a study of Ca n emission lines from SN 
1987 A, Li & McCray (1993) recognized that this feature must 
result from fluorescence following absorption of UV radiation 
by the Ca ii AA3934, 3969 (H, K) lines. This result is confirmed 
by our Monte Carlo calculations. 

The dashed vertical marks in the top rows of Figures 4 and 5 
indicate some of the most prominent of the thermally excited 
forbidden lines of Fe ii, and these are also listed in Tables 1 and 
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Fig. 5. — The emergent spectrum for a UHeS model at 800 days and 2000 K 


TABLE l 

Strong Emission Features on 
Day 300 at 6000 K 
A. Fluorescence Lines 


Elements and Their 
A (A) Contributions (%) 


3765 Ti ii( 33), V n(33), Fe n(22) 

4037 Ti n(52), V u(40) 

4640 Fe n(45), Ti ti(32), Cr ji(23) 

6254 Sc n(49), Fe n(44) 

6508 Fe n((57), Ti n(35) 

8514 Ca ii(100) 

8685 Ca 11(100) 

10493 Fe ii(90), Ti ii(8) 

10811 Fe ii(99) 

30733 Feu(100) 


B. Strong Fe a Forbidden Lines (A or /an) 


3377 

3506 

4289 

4415 

4729 

5275 

7455 

2.05 


3442 

3540 

4321 

4418 

4816 

5335 

8619 

2.13 


3453 

4178 

4348 

4453 

4891 

5378 

8894 

17.94 


3456 

4245 

4354 

4459 

4907 

7157 

1.26 

24.52 


3503 

4278 

4360 

4476 

5160 

7174 

1.64 

25.99 


3505 

4281 

4361 

4490 

5263 

7390 

2.01 


TABLE 2 

Strong Emission Features on 
Day 800 at 2000 K 
A. Fluorescence Lines 


Elements and Their 
A (A) Contributions {%) 


2766 Fe n((87), Cr n(10) 

2965 Fe n(63), V n(36) 

3276 Ti ii(71), V u(20), Fe ii(10) 

3765 Ti ii(93) 

4328 Sc n(64), Ti ii(32) 

4415 Sc ii(51), Ti n(42) 

4926 Fe n(92) 

5025 Fe u(84) 

5177 Fe ii(87) 

8514 Ca 11(100) 

8685 Ca n(89), Mn u(ll) 

30733 Fen(lOO) 


B. Strong Fe n Forbidden Lines (A or ^m) 


7157 

7174 

7390 

7455 

8619 

8894 

1.26 

1.29 

1.53 

4.08 

4.12 

4.61 

4.89 

5.34 

17.94 

25.99 




rrnr 


UV OPACITY AND FLUORESCENCE 


377 


TABLE 3 

Energy Budget (ergs s“ ') 


Luminosity 

300d, 6000 K 

800d. 2000 K 

L(initial photons) 

6.3 x 10 38 

6.1 x 10 36 

^emergent photons) 

3.6 x 10 38 

4.2 x I0 3h 

L(emergent)/L(initial) (%) — 

56 

68 

L(Fe ii forbidden lines) 

2.9 x 10 39 

7.9 x 10 36 


2. In principle, it should be possible to distinguish spectro- 
scopically the contributions of thermal emission and fluores- 
cence to the optical radiation from the supernova envelope. 

Table 3 lists the energy budgets for the models represented 
by Figures 4 and 5. The first row gives (for each model) the 
assumed luminosity in the He i two-photon continuum (for 
k > 912 A), and the second row gives the luminosity of the 
emergent fluorescence spectrum, not including forbidden line 
emission. The emergent luminosity is less than the input lumi- 
nosity because a significant fraction (see Table 3, third row) of 
the input luminosity pumps the populations of metastable 
levels. However, it is better not to regard the forbidden lines 
resulting from the decay of these metastable levels as fluores- 
cence because the populations of these levels are largely con- 
trolled by collisions with thermal electrons, not fluorescent 
pumping. Indeed, the total luminosity of the Fe n forbidden 
lines, listed in the fourth row, exceeds the difference between 
the luminosity of the input two-photon continuum and the 


emergent fluorescent spectrum by a factor 10.7 at 300 days, and 
by a factor 4.1 at 800 days. Therefore, UV pumping of the 
metastable levels cannot heat the gas enough to account for the 
forbidden line emission. This is not a problem for the model, 
however, because fast electrons will deposit a greater fraction 
of their energy as direct heating of the electron gas than as 
excitation of two-photon continua (Xu & McCray 1991; 
Kozma & Fransson 1992). 

We also simulated a few other models to test the sensitivity 
of the emergent spectrum to the assumptions of the model. For 
example, Figure 6 shows the emergent spectrum for a model 
with an envelope identical to that of Figure 5 (r — 300 days, 
T = 6000 K), and a He i two-photon continuum with the same 
intensity, except that the source is now a central point source 
instead of a uniformly distributed source. As one would expect, 
in Figure 6 the UV spectrum is more strongly absorbed and 
the optical fluorescence is stronger than in Figure 4 but other- 
wise the emergent spectra look very similar. 

Figures la and lb may provide some insight into this behav- 
ior. Here, each point represents one emergent photon selected 
from a single Monte Carlo trial. The abscissa indicates the 
wavelength of the original UV continuum photon, and the 
ordinate indicates the wavelength of the emergent photon. In 
each panel the points lying on the diagonal line represent con- 
tinuum photons that escape without scattering. Note that the 
diagonal line is well-covered with points in the uniformly dis- 
tributed source model (panel a ), but that very few photons with 
k < 3000 A escape freely in the central source model (panel b). 



0.1 0.2 O.a 0.4 0.5 0.7 1 2 0 4 

Wavelength (/Am) 

Fig. 6. — The emergent spectrum for a central He i two-photon continuum source (CHeS) at 300 days and 6000 K 
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Fig. 7.— Wavelengths of the emergent photons vs. those of the initial photons. Panels {a) and ( b ) are for the cases of Fig. 4 and 6, respectively. 


Evidently, the UV photons that escape freely are those created 
near the surface. The fact that the emergent photons cluster in 
horizontal bands indicates that the emergent fluorescent spec- 
trum is rather insensitive to the input spectrum — its shape is 
determined more by atomic branching ratios. The points 
below the diagonal line are examples of events in which a 
continuum photon is absorbed by a resonance transition from 
a thermally excited metastable level, after which the excited 
atom decays to a lower metastable level. 

To gain further insight into the scattering and fluorescence, 
we show in Figures 8 a and 8 b the integrated probability dis- 
tributions, P(n\ that an emergent photon will escape after n or 
fewer scatterings. (These distributions are derived from the 
results of step [7] in § 3.) The solid curves represent models 
with uniformly distributed source photons, whereas the dashed 
curves represent models with central sources. In Figure 8a the 
envelope is assumed to have a temperature 6000 K, and the 
dashed and solid curves correspond to the models represented 
by Figures 4 and 6, respectively. Note that approximately 20% 
of the initial photons escape without scattering, whereas 80% 
of the photons escape with 1 7 or fewer scatterings in the dis- 
tributed source model, and with 25 or fewer scatterings in the 
central source model. 

The effect of reducing envelope temperature to 2000 K is 
illustrated by Figure 8 b. The curves resemble those in Figure 
8a, except that the characteristic number of scatterings 
required before a photon escapes is reduced by a factor ~ 5, 
owing to the greater number of transparent windows (compare 
Figs. 3a, b). 


Figures 9 a and 9b show the integral probability distributions 
for the emergent photons as a function of last scattering radius. 
With a uniformly distributed source, the shapes of the prob- 
ability distributions ( solid curves) represent in first approx- 
imation the r 3 dependence of the integrated volume. More 
interesting are the dashed curves, showing the probability dis- 
tributions for central source models. They show that there is 
considerable propagation in radius before a typical photon 
finally escapes. Thus, the escape process is a combination of 
migration in wavelength and spatial diffusion. 

In our calculations we have neglected overlap of resonance 
lines. In fact, however, the resonance shells have finite thick- 
ness, given by equation (1). If the resonance shells of two ore 
more lines ovelap, and one or more of the overlapping lines has 
substantial optical depth, our approximation that the reso- 
nance shells are infinitely thin will introduce an error in the 
resulting fluorescence spectrum. Our approximation will 
convert too much of the absorbed UV radiation into the fluo- 
rescence spectrum of the bluer component at the expense of the 
fluorescence spectrum of the redder component. The approx- 
imation will not significantly affect the effective UV opacity, 
however, because the net blocking due to two overlapping lines 
is the same as that due to two nearby separated lines. 

To see whether overlapping lines might introduce a signifi- 
cant error, we counted the number of resonance lines with 
Sobolev optical depth above a given threshold, and also the 
number of pairs of such lines with wavelengths separated by 
less than the thermal width. We find, for example, that at 300 
days and 6000 K, 6067 Fe n lines have optical depth > 
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Fig. 8 —Fractional number of the photons vs. number of scatterings, n scl . (a): day 300, 6000 K ;(b): day 300, 2000 K. Solid: UHeS mode!; dashed: CHeS model. 


0.1, and that only about 3% of these lines overlap. Therefore, 
we conclude that our neglect of line overlap will introduce 
errors of at most a few percent in the fluorescent spectrum. 


6. DISCUSSION 

Here, we first comment on the relationship of the present 
study to other studies of resonance line transfer and then 
discuss future work to be done. 

Beginning with the pioneering work by Karp et al. (1977), 
many authors have studied the opacity of expanding super- 
nova envelopes due to multiple line scattering. Most of these 
studies have been directed primarily at understanding super- 
nova spectra during their early photospheric phases. Wagoner 
et al. (1991) presented an improved formulation of the effective 
continuum opacity due to many resonance lines, using an 
extensive database of atomic transitions compiled by Kurucz. 
Wagoner et al. assumed that both the ionization balance and 
the level populations were in LTE, and modeled resonance line 
emission and absorption with a two-level atom approximation, 
in which true absorption and emission were characterized by a 
thermalization parameter. Lucy (1987) showed that the UV 
radiation from SN 1987 A at early times was strongly blocked 
due to resonance scattering by many lines of metals, a result 
that was confirmed in more detailed model atmosphere calcu- 
lations by Eastman & Kirshner (1989), Schmutz et al. (1990), 
and Hauschildt & Ensman (1994). None of these calculations 
has attempted to model the effects of fluorescence in the trans- 
fer of UV radiation. For example, Hauschildt & Ensman (1994) 


parameterized the UV line scattering with an effective albedo 
that is independent of wavelength. 

As we have seen, during the nebular phase of SN 1987A the 
ionization states and level populations of metals are far from 
LTE. UV radiation created by gamma-ray energy deposition 
in a supernova envelope will be largely blocked by resonance 
scattering in the envelope. Roughly 60% of the blocked UV 
luminosity will be converted by fluorescence into a quasi con- 
tinuum of thousands of optical and infrared emission lines. 
This quasi continuum resembles the continuum seen in the 
optical spectrum of SN 1987 A at late times. Moreover, we have 
shown that the UV opacity of an expanding supernova 
envelope will decrease with time as a result of the decreasing 
temperature and density of the envelope. This provides a 
natural explanation for the fact that the fraction of the bolo- 
metric luminosity of SN 1987 A that emerges as UV radiation 
increases with time (McCray 1993; Pun et al. 1995). The model 
calculations that we have presented here give us confidence 
that we are well on the way to providing a quantitative and 
detailed account for the time evolution of the entire spectrum 
of SN 1987A. 

We hasten to add, however, that the calculations presented 
here are only a step on the way to building a model for the 
spectrum of SN 1987A. To reach that goal, we must consider a 
more realistic model envelope for SN 1987A, including physi- 
cal processes not considered here. The tasks to be completed 
are as follows. 

1. We must consider a model envelope with a more realistic 
distribution of composition and density. In the present models, 
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Fig. 9. — Fractional number of emergent photons vs. radius of last scattering, (a): day 300, 6000 K; (6): day 300, 2000 K. Solid: UHeS model; dashed: CHeS 
model. 


we have assumed a uniform density hydrogen envelope with 
cosmic abundance ratios of trace elements. In fact, at early 
times a substantial fraction of the gamma-ray luminosity will 
be deposited in the newly synthesized Fe/Co/Ni bubbles (Li et 
al. 1993), so we must consider the production and fluorescence 
of UV photons in these bubbles. Likewise, we expect that com- 
parable fractions of the gamma-ray luminosity are deposited in 
the hydrogen-rich and helium-rich regions of the inner 
envelope, and so we must consider these regions separately. 

2. We must include the contribution of other species, not 
considered here, to the UV opacity and fluorescence. For 
example, at late times, and in the Fe/Co/Ni bubbles, we cannot 
be sure that all the metals are once-ionized. Therefore, we must 
consider the ionization balance of metals more carefully and 
include the possible contributions of neutral metals. Likewise 
we must consider the contributions of molecules to UV opacity 
and fluorescence. For example, Culhane & McCray (1995) 
have shown that about 1% of the hydrogen in the inner 
envelope will be molecular, and that the Lyman and Werner 
bands of the H 2 molecules will certainly make a major contri- 
bution to the UV opacity and fluorescence. 

3. We must include more realistic models of the UV source 
spectrum. In the models presented here, we only illustrated the 
degradation and fluorescence of the He i two-photon contin- 
uum. But gamma-ray energy deposition will also produce a H i 
two-photon continuum and a spectrum of UV emission lines 
excited by the stopping of nonthermal electrons. 

4. Since we have shown that the UV opacity and fluores- 
cence are very sensitive to temperature, we must model the 
evolution of temperature in the envelope more carefully. 


5. In making fits of model spectra to observations, we must 
include other sources of continuum opacity. For example, Xu 
et al. (1992) showed that Balmer continuum opacity in the 
hydrogen envelope would be substantial for t < 200 days. 
Moreover, extinction of optical and UV radiation by dust 
grains in the envelope is known to be substantial for t > 400 
days. 

This list is substantial, but we are making steady progress 
toward completing these tasks. We think the effort is worth- 
while for a number of reasons. First, our attempts to model the 
evolution of the spectrum of SN 1987 A often point the way to 
new spectral diagnostics of the physical conditions in the 
envelope. (An example from the present study is the ability to 
distinguish contributions to the optical spectrum from thermal 
emission of forbidden lines and by fluorescence, which may 
provide a sensitive measure of temperature.) Second, SN 
1987A is truly a Rosetta stone for supernova spectra. With a 
detailed understanding of the physics involved in the forma- 
tion of various features in the spectrum of SN 1987 A, we can 
more confidently draw physical inferences from such features 
in the spectra of other supernova, for which such com- 
prehensive observations and physical models are lacking. 
Third, the models elucidate properties of resonance line trans- 
fer and fluorescence that are important not only for super- 
novae, but also for classical novae and stellar winds. 


This work was supported by NASA through grant NAGW- 
2900. 
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AN ANALYTIC TEST 


To test our Monte Carlo code we consider the following idealized model. Assume that a uniform, isotropic stellar envelope of 
radius R is expanding homologously with a maximum expansion velocity (at R), u max = 0.01c. Assume that there is only one 
artificial element in the envelope. This element has several hundred equally spaced excited levels, each of which has only one 
resonance transition, to the ground level (no subordinate transitions and no fluorescence). The wavelengths of these lines are given 
by 

A, = / 0 + iSX , i = 1,2,3..., (Al) 

where SX is given by 


<5/l = 0.U o — . (A2) 

c 

Assume that all of these lines are optically thick, so that any photon that comes into resonance with one of these lines must be 
scattered. Therefore, if a photon with original wavelength X 0 appears at the center of the envelope, it must be scattered at least 10 
times before it reaches the surface, and it should emerge with wavelength X > X 0 {\ + v m Jc). The minimum number of scatterings 
corresponds to the case where every scattering is in the forward direction. 

We assume, however, that the scattering angle is random. Since the photon’s wavelength must increase by 6X after each scattering, 
we can write the wavelength of the emergent photons as a function of the number of scatterings, n scl 

*("sct) = ^0 + "set ^ • < A3 ) 

The propagation of the UV photons should be described accurately by the diffusion equation 

= D V 2 N(r, t) - P abs N(r, f) , (A4) 

dt 

where N(r , t ) is the total number density of photons at distance r from the center, t is the diffusion time in units of mean free travel 
time (i.e., t = n scl ), P abs is the probability of absorption per scattering, and D = 5 is the diffusion coefficient. (There is a subtle point 
here. For diffusion through a random distribution of targets, D = 3, because <(<5r) 2 > = 2 (Sr > 2 , where Sr is the distance traveled per 
scattering. But in the case at hand, ((Sr) 2 ) = (Sr) 2 .) 

Consider the problem of a central source of UV photons. The boundary conditions are 

V(0, t > 0) (A5) 


is finite, and 

N(R\ t) = 0 . 

Following Illarionov et al. ( 1 979), we set R' = R + R/(3/2 R - 1 ) in equation (A6). The initial condition is 

N(r, 0) = -^ S(r ) , 

4nr~ 


(A6) 


(A7) 


where S(r) is the Dirac delta function, and Jf is the total number of photons emitted at the center. Following Illarionov et al. (1979), 
we define a new variable /(r, t) = rN(r , f), so that equation (A4) becomes 


df(r, t) d 2 f(r, t) 

dt dr 2 


~ Pahs /(>% 0 • 


(A8) 


The general solution to equation (A8) is 

/ <r, t) = ~ exp ( - P abs t) X n exp 
IK 1 



(A9) 


and the luminosity emerging at the true surface at time t is 

An 

L(R , r) = — 


_/('• '» “ ' Tr /(r ' 


(A10) 


This luminosity is equivalent to the emergent spectrum because the wavelength of the emergent photons is determined by the 
number of scatterings (eq. [A3]), which is proportional to the travel time of the photons within the envelope. Therefore, we can write 
the emergent spectrum as 

DC 2 00 

F(A) = .A r 2n — y exp ( - P abs n sct ) X n ex P 

K n = 1 


nVC 2 

R 2 


Dn scl )[sin (miQ — nn( cos (nn()'] 


(All) 
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Fig. 10. Spectrum of emergent photons (top) and absorbed photons (bottom) for an artificial element in a uniform supernova envelope. Solid curves • analytic 
results; diamonds : results from Monte Carlo simulation. 


where 


, (A12) 

C = R/R’- (A 13) 

We can also calculate the flux of the absorbed photons. The survival probability per scattering is given by P sur (l) = j _ p abj 
Thus, after n scatterings, the survival probability is P s Jn) = (1 - P abs )" and the absorption probability at the [(n +"i)th] scattering 
is given by P abs (n + 1) = P sur (»)P abs (l) = (1 - P ab J"P al „. P abs (n) is in fact the measure of the flux of the absorbed photons. So, we can 
write this flux 

^absW = (l“^b 8 r C ‘ a) / > abs, (A 14) 

where n sct is a function of k given by equation (A 12). 

To test our Monte Carlo code, we consider 10 5 initial photons with wavelength k 0 = 1000 A emitted at the center of the envelope, 
and assume P abs = 0.05. The resulting spectrum of the emergent and absorbed photons are shown in Figure 10, top and bottom’ 
respectively. The diamonds are the Monte Carlo results and the solid curves are the analytic curves from equations (All) and (A 14), 
respectively, The excellent agreement confirms the diffusion model presented above and the reliability of our code. 
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